A "Burnt Bridge" Brownian Ratchet 
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Motivated by a biased diffusion of molecular motors with the bias dependent on the state of 
the substrate, we investigate a random walk on a one-dimensional lattice that contains weak links 
(called "bridges") which are affected by the walker. Namely, a bridge is destroyed with probability 
p when the walker crosses it; the walker is not allowed to cross it again and this leads to a directed 
motion. The velocity of the walker is determined analytically for equidistant bridges. The special 
case of p = 1 is more tractable — both the velocity and the diffusion constant are calculated for 
uncorrelated locations of bridges, including periodic and random distributions. 

PACS numbers: 02.50.-r, 05.40.Fb, 87.16.Nn 



INTRODUCTION 



The motion of a particle depends on the medium. Of- 
ten the inverse is also true, that is, the particle motion 
changes the medium. Such problems arc characterized 
by an infinite memory — not only the present position 
of the particle, but the entire past determines the future 
— and they are usually extremely difficult. Perhaps the 
most famous example is the self-avoiding walk which is a 
random walk on a lattice with the restriction that hops 
to already visited sites are forbidden yj. Similarly in a 
path-avoiding walk, a random walker is not allowed to go 
over already visited links. Here we examine a generaliza- 
tion of the path-avoiding walk in which the medium is a 
lattice with two kinds of links, strong and weak: strong 
links are unaffected by the walker while weak links, called 
bridges, "burn" when they are crossed by the walker. The 
random walker is not allowed to cross a burnt bridge. 
Obviously the burnt bridge model reduces to the path- 
avoiding walk if all links are weak. We also investigate 
a stochastic burnt bridge model in which the crossing 
of an intact bridge leads to burning only with a certain 
probability p, while with probability 1 — p the bridge re- 
mains intact. 

The stochastic burnt bridge model is a simplification 
of models proposed to mimic classical molecular motors 
[J 13 0] with energy coming from ATP hydrolysis • 
Recently it has been shown [8| that the stochastic burnt 
bridge model with two weakly coupled tracks (the walker 
moves on the ladder and the hopping rate between the 
tracks is small compared to the hopping rate along the 
tracks) accurately describes experimental results on the 
motion of an activated collagenase on the collagen fibril. 

We shall focus on the one track stochastic burnt bridge 
model. Forbidding the crossing of burnt bridges essen- 
tially imposes a bias, and the goal is to compute the ve- 
locity v(c,p) and the diffusion coefficient D(c,p) as func- 
tions of the density of bridges c and the bridge burning 



probability p. The velocity and the diffusion coefficient 
also depend on the positioning of the bridges. We shall 
tacitly assume that the bridges are placed without corre- 
lations, and we shall often specify our findings to two 
particularly interesting and natural positioning of the 
bridges — a regular equidistant spacing and a random 
distribution. 

The rest of this paper is organized as follows. In the 
next section, we describe various versions of the burnt 
bridge model and outline the major results. Section ITTT1 
is devoted to the derivation of the velocity and the diffu- 
sion coefficient for the burnt bridge model and a modi- 
fied burnt bridge model. In Sec. lIVI the stochastic burnt 
bridge model is studied, and the velocity is computed for 
equidistant bridges. Finally, a few open questions are 
discussed (Sec. EJ. Various calculations are relegated to 
the Appendices. 



II. MAIN RESULTS 

In all versions of a burnt bridge model, the walker 
asymptotically behaves as a biased random walk. Math- 
ematically this means that the probability of finding 
the walker at position a; is a Gaussian centered around 
(x) = vt with width (x 2 ) - (x) 2 = 2Dt. 

Away from the bridges, the walker hops to adjacent 
sites equiprobably (Fig. [I}. Thus if there were no bridges, 
the velocity would be equal to zero and the diffusion con- 
stant would be equal to 1/2 (the lattice spacing and the 
time step between the hops are set to unity). Bridges 
(which are assumed to be intact initially) generate a di- 
rected motion. The velocity v of the walker depends on 
the density of bridges c and on the bridge distribution. 
Particularly simple expressions are obtained for periodi- 
cally and randomly positioned bridges 



v{c) 



c periodic 
c/(2 — c) random 



(1) 
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Of course, in the periodic case the density c attains only 
inverse integer values (1, 1/2, 1/3, . . .) while when bridges 
are placed at random the density can attain any value 
< c < 1. 
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In the realm of continuum approximation, the veloc- 
ity v(c) has been computed by Mai et al. @- Unex- 
pectedly, the continuum approximation is exact in the 
periodic case; for random locations, the continuum ap- 
proximation gives v(c) — c/2. This result is asymptoti- 
cally exact in the small c limit, albeit overall it is just an 
approximation which more and more deviates from the 
exact result as the density approaches c = 1 (that is, for 
the path-avoiding walk) . 




FIG. 1: The random walker (filled triangle) hops to adjacent 
sites when it is away from the bridges. Strong links are shown 
by horizontal lines, intact bridges are depicted by arcs, and 
the absence of a link indicates a burnt bridge. 




FIG. 2: Velocity of the walker in the original (solid) and in 
the modified (dashed) burnt bridge model, for random and 
regular bridge distributions. 



The dependence of the diffusion coefficient on the den- 
sity of bridges c is also relatively simple for periodically 
and randomly positioned bridges 
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(1 



D(c) 



2(1- c/2) 2 



periodic 



random 



(2) 



The diffusion coefficient monotonously decreases as c in- 
creases (see Fig.|3J), and D(l) = since on a lattice fully 
covered by bridges, the walker moves deterministically. 
The diminishing of D{c) has apparently been observed in 
Ref. [3. Intriguingly, equation (J2J gives D per (+0) = 1/3 
(-Dran(+0) = 3/2) which is smaller (larger) than the 
"bare" diffusion coefficient Dbarc = 1/2 that charac- 
terizes diffusion on the one-dimensional lattice without 
bridges. This sudden jump of the diffusion coefficient 
occurs when the density becomes positive. The reason 
is that any positive c (irrespectively however small it is) 
makes the lasting influence on the fate of the walker who 
is forced to remain to the right of the last burnt bridge. 
Thus the very rare burning events substantially affect the 
diffusion coefficient. 

The details of the walker dynamics at the boundary 
of the burnt bridge affect the velocity and the diffusion 
coefficient. To illustrate this assertion recall that in the 
framework of the burnt bridge model the walker at the 
boundary of the burnt bridge always moves to the right. 
Another natural definition is to allow an attempt to cross 
the burnt bridge — the attempt fails and the walker re- 
mains at its position. We calculated the velocity for this 
modified burnt bridge model 9]: 



v(c) 




c) periodic 
random 



(3) 




FIG. 3: Diffusion coefficient in the original (solid) and in 
the modified (dashed) burnt bridge model, for random and 
regular bridge distributions. 



We also computed the diffusion coefficient: 
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— — periodic 
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Perhaps the largest difference between the two models is 
that in the realm of the modified burnt bridge model the 
walker never moves deterministically — even when c = 1 
it moves diffusively although the diffusion coefficient is 
small, namely it is 4 times smaller than the bare diffusion 
coefficient. Not surprisingly, the quantitative predictions 
of the two models are substantially different when c is 
large (see Figs. |21 and I3J). 

For the stochastic burnt bridge model, we computed 
only the velocity, and only for periodically located 
bridges. We found 



v(c,p) 



cp 



2-p+U 

cp + 2-p p(l-c) + V 



(5) 
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FIG. 4: The velocity v(c,p) vs. c (for p = 1, 1/2, 1/4, 1/8, top 
to bottom) and p (for c = 1,1/2,1/4,1/8, top to bottom). 
The lattice is periodically covered by bridges. The symbols 
are the results of numerical simulations. 
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FIG. 5: The fraction of intact bridges, I(c,p) vs. c (for p — 
1/4,1/16,1/64,1/256, from bottom to top), and vs. p (for 
p — 1, 1/4, 1/16, 1/64, from top to bottom) in the periodic 
case. 



where we used the shorthand notation 
p(2 -}»)(! -c) 



V 




Ac 



p(2-p)(l-c) s 



For p = 1, Eq. JSJ agrees with already known result 
v = c [sec Eq. Q]; for c = 1 (the lattice fully covered 
by bridges), the velocity is given by the following neat 
expression 



v(l,p) 



p + x /p(2 - p) 



(6) 



From Eq. (J5J (see also Fig. one finds the asymptotics 



v(c,p) 



when c«p 



I y/ cp/2 when p«c 



(7) 



For c <§C p, the distance between neighboring bridges is 
large. Thus the walker typically crosses the next bridge 
several times and hence almost all bridges get burnt. 
Therefore the p = 1 results ought to be recovered. Equa- 
tion Q shows that this is indeed correct in the periodic 
case; in the random case (where we do not know an exact 
solution) we similarly expect v(c,p) m c/2 when c <C p. 
In the complimentary limit p <C c, the walker on average 
makes many steps before the burning occurs, and it is 
intuitively obvious that we can renormalize c — > 1 and 
simultaneously p — > cp. Hence v(c,p) —> u(l,cp) when 
p <C c, and therefore the asymptotics given in JjJ can 
also be extracted from the simple solution (jfjj). 



Another interesting and experimentally accessible 
quantity is the fraction of bridges left intact by the 
walker. In the long time limit it can be expressed as 



I(c,p) = l-V- CP + 2 - p 
V ,FJ c V + 2-p 



(8) 



Figure |3] shows that the fraction of intact bridges is a 
decreasing function of p for fixed c (this is intuitively 
obvious) and an increasing function of c for fixed p. This 
latter feature is understood by noting that on average a 
bridge is visited more often when the density of bridges 
gets smaller. 



III. THE BURNT BRIDGE MODEL 

The walker hops x — > x ± 1 equiprobably when it is 
away from the burnt bridges. The lattice spacing and the 
time step between successive hops are set to unity, and 
therefore the diffusion coefficient is D — 1/2 for the lat- 
tice with no bridges. We assume that initially all bridges 
are intact. Let the walker cross a bridge for the first time 
from the left. This implies that the walker will never cross 
it again and determines the fate of the walker, namely the 
walker will drift to the right and the closest burnt bridge 
will always be on the left. By definition, when the walker 
is on the right boundary of the burnt bridge, it always 
makes the step to the right. 

Two special positioning of the bridges are regular 
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(equidistant locations) and random. These are two most 
interesting cases, although the formalism can be straight- 
forwardly extended to arbitrary distributions as long as 
the distances between neighboring bridges are uncorre- 
cted. 



A. Velocity and Diffusion coefficient 

The walker successively crosses the intact bridges 
which immediately turn into burnt bridges. Instead of 
the original walker it proves convenient to consider an 
equivalent walk whose position at time t is on the right 
edge of the the last burnt bridge. This walk is still dis- 
crete in space and time, but the step length is now equal 
to the distance between successive bridges, and the time 
between the steps is at least as large as the step length. 
In general, we denote by &(x, t) the probability that the 
walk makes a step of length x after waiting t units of 
time. 

The following derivation is essentially the discrete time 
version of the so-called continuous time random walk 
(CTRW) [l(J El This title refers to the continuous na- 
ture of the the waiting time distribution, while in our 
model the discrete nature of time is essential, and there- 
fore we present the complete derivation of the necessary 
results. Note also that even though our walk always steps 
to the right, the derivation is the same as for a general 
walk which can step in both directions; hence, we present 
a general derivation. 

Let ^>j(x,t) be the probability that the walk arrives at 
site x at time t and at the j'-th step, and hence (x, t) = 
^(x,t). This probability satisfies the recurrence formula 



9j(x,t)= J2 X>;-i(s *')*(*>*) • (9) 

x' — — oo t'—O 

Using the generating function (the discrete version of the 
Fourier-Laplace transform) 



oo oo 



V (10) 



x — — oo t— 

we transform the convolution @ into the product 

9j(q,u) = *j_i(q,u)*(q,u) = [V(q,u)} j . (11) 

The probability Q(x,t) of arriving at site x at time t 
(after an arbitrary number of steps) is 



(12) 



The corresponding generating function can be written in 
a closed form 



oo 1 



3=0 



(13) 



The probability P(x,t) that the walk is at site x at 
time t can be obtained by noting that in order to be at 
site x at time t, the walk has to arrive at site x not later 
than at time t and has to stay there until t. Hence 

t 

P(ar,t) = 530(t-OQ(»»* / ) (14) 

t'=0 

where <p(t) is the probability that the walk does not move 
for a time interval t 



4>{t) 



(15) 



The probability ^(t) that the walk makes at least one 
step during the time interval t is obtained by summing 
over all possible step lengths 



(16) 



Using Eq. I|15|l we compute the generating function of 



<j>(u) 



1 - #(u) 



(17) 



The generating function of P(x, t) is now easily derived 
since Eq. (|14|) is a convolution: 



P{q,u) 



(u)Q(q,u) = 



1 - 



(l-u)[l -*(«,«)] 



(18) 



We are going show that position of the walk follows a 
Gaussian probability distribution for large times 



P(x,t) 



1 



exp 



{x - vt) 2 
ADt 



(19) 



which is centered around a mean value (x) = vt with a 



mean square deviation (x 



2Dt, with v being 



the speed of the walk, and D being the diffusion coeffi- 
cient. We define the new variables 7 and e as q = e 



?7 



and u = e~ e . In order to show that the long time limit of 
P(x, t) is Gaussian, we ought to show that in the 7, e — > 
limit P(7, e) is equal to the Laplace-Fourier transform of 
the Gaussian 119fl 



P(7,e) = 



ijv + r y 2 D 



(20) 



Therefore we must show that P(7, e) given by Eq. i|18|) 
attains the form of Eq. (|20|l . 

Although the probability \&(x, i) is not separable in 
general, it can always be written as the product 



®(x,t) = S(x)V(t\x) 



(21) 



of the probability S(x) that the next step has length x 
(distance to the next bridge) times the conditional prob- 
ability fy(t\x) that the next step happens after t waiting 
time, given that the length of this step is x. 
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Now we calculate ^(7, e) up to the second order in 7 
and e. First, we calculate the generating function with 
respect to time in the e — > limit. Plugging H21(l into 

*(x,e) = ^*(M)e- £ * 



and expanding in e up to the second order we obtain 



where the moments of time are calculated at some fixed 
x length of interval 



t=0 



(22) 



Performing the Fourier transform of ^{x,u) and taking 
the 7^0 limit, we arrive at 



*( 7) c) = l-e{[t] x )+i 1 {x)^e 1 (x[t] x ) + -{[t%)— '-(x 2 ) 

up to second order in 7 and e. We also need to calculate 
^(u), the generating function of *&(t) defined by Eq. (|16[l . 
It is sufficient to know it only up to first order in e: 

*(e) = 1 - e([t] x ) 

Plugging the above results into Eq. I |18|l we find that 
up to first order in both 7 and e, the quantity \I/(7,e) 
attains the form P(7, e) = (e — ijv)~ x with 



(23) 



Thus P{x, t) is centered around (x) — vt with velocity 
given by Eq. 

Calculating P("f, e) up to second order, and using the 
first order expression for e w —ijv in the terms contain- 
ing £7 and e 2 , we arrive at Eq. I|2U|) with the diffusion 
coefficient given by 



D 



2<[*] a 



([t%)(x) 2 (x[t] x )(x) 

2<w,) 3 (m 2 



(24) 



Thus we conclude that P(x, t) is indeed Gaussian. Note 
that this result applies to any random walk — discrete or 
continuous — where the steps are uncorrelated and all of 
the moments used in Eq. (|24|l are finite. Specifically, it 
applies to the burnt bridge model (p = 1) if the distances 
between bridges are uncorrelated. 



1. Special cases 

Consider first equidistant bridges separated by dis- 
tance £. Then S(x) = S x j, and therefore (x) = £, (x 2 ) — 



£ 2 , etc. The velocity |j53J| and diffusion coefficient J2IJ| 
simplify to 



[t]t ' 2 



m 

Using the moments of time computed in Appendix lAl 
[Eqs. (|A5|1 and (|A9|) ]. we arrive at Eqs. and 10). 

For randomly distributed bridges, the probability that 
the walk makes a step of length x, that is the probability 
of having two neighboring bridges at distance x > 0, is 



S{x) = c(l - cf 



(25) 



We again use Eqs. (|A5|) and l|A9(l for the moments of 
time, and we also need the first four moments of x 



1 



[x) = - 
c 



2-c 



6 — 6c + c 2 



(26) 



24 - 36c + 14c 2 - c 3 



Using these expressions in Eqs. l{2"3"|) and (J52JI, we obtain 
the velocity and the diffusion coefficient given by Eqs. 
and J2J), respectively. 

Finally, consider the bimodal distribution 



S(x) = qi5 x ^ + q 2 6 X ,£ 2 



(27) 



with two possible separations between the bridges, £\ and 
£2, occurring independently with respective probabilities 
91 and 92 (of course, q±, q 2 > and q\ + q% = 1). In this 
situation, the velocity l|23|) becomes 

_ 9l^l + 92^2 
V ~ qi£\ + qA 

and the diffusion coefficient (|24() turns into 



D = 



1 _ (gl^ + <Z2^)(<Zl4+<^ 2 ) 

2 
5 
6 



( qi £ 2 + q 2 £ 2 ) 2 
q2£t){qi£l + q2£2? 



1 (gill 

3 M 2 



92^2 



6(91^+92^) 



2\3 



92^) 2 
(28) 



For the system with bimodal bridge distribution l|27|l . 
and apparently for an arbitrary uncorrelated positioning 
of bridges, the diffusion coefficient exceeds that of the 
corresponding periodic system at the same bridge den- 
sity. This general assertion is easy to verify in a partic- 
ularly interesting case when the density of bridges van- 
ishes. Taking the limit £\, £ 2 — > 00, and keeping the ratio 
£i/£2 = r(< 1) constant, we recast Eq. <|28|) into 



D= X - 
2 



/ qir-\ 


- q2 x 


2 

) 


"5 


4 

q\r 


+ 92 


qi r 3 + g 2 


\qir 2 - 


f 92, 




6 


qir 2 


+ 92 


qir + 92 



(29) 
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A straightforward analysis of Eq. I|29[l shows that the dif- 
fusion coefficient is larger than 1/3, which is the diffusion 
coefficient in the periodic case (qi = or q2 — 0). From 
Eq. I|29l) one finds that for a -C 1, the maximal diffu- 
sion coefficient, approximately D ss gf r_2 ' 1S achieved 
when the system is predominantly composed of shorter 
segments l\, namely when q-i m r 2 /2. Therefore a "su- 
perposition" of two equidistant distributions, each char- 
acterized by the diffusion coefficient 1/3, may have an 
arbitrarily large diffusion coefficient. 



2. Simulations 

The velocity and the diffusion coefficient of the walker 
are determined using the basic formulas v = (x) jt and 
D = ((x 2 ) — (x) 2 )/2t. Since the motion is self-averaging, 
simulating a single walker for a long time is sufficient to 
obtain the velocity. To measure the diffusion coefficient, 
however, one has to perform averages over several runs. 
In the case of randomly distributed bridges, one also has 
to average over the bridge distribution. 

Figure shows numerical results for the diffusion co- 
efficient at various times. The convergence to the the- 
oretical predictions is slow when the density of bridges 
is small. During a short time interval the walker does 
not reach the second bridge, and actually behaves as a 
simple random walk with a reflecting boundary at the 
origin. Hence the probability of finding the particle at 
position x > is a Gaussian centered around the origin, 
and the formal definition of the diffusion coefficient yields 
D = 1/2 — 1/-7T. For the time intervals large compared to 
the time (of the order of c -2 ) between overtaking succes- 
sive bridges, the coarse-grained motion becomes similar 
to a biased random walk with the diffusion coefficient 
approaching the theoretical predictions: D(+0) = 1/3 in 
the periodic case and D(+0) — 3/2 in the random case. 
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FIG. 6: The diffusion coefficient D(c) for the random and 
periodic bridge locations. Results of the simulations are also 
displayed for comparison. The arrow points to the theoretical 
value D = 1/2 — 1/tt corresponding to a random walk with a 
reflecting boundary. 



B. Correlation function 

A correlation function measured experimentally in 
Ref. [j| is apparently proportional ^3 to the probability 
C(t) that the walker at the site xq, will be at the same 
position time t later. 

As a warm-up, consider the extreme cases of the lat- 
tice without bridges (c = 0) and the lattice fully covered 
by bridges (c = 1). In the former case, the correlation 
function obviously vanishes for odd t while for even t it 
is given by the well-known expression 

'2t> 



C(2t) 



-21 



(30) 



Note that the correlation function decays algebraically in 
the large time limit: 

C(2t) -> as t -> oo (31) 



7rf 

For the lattice fully covered by bridges the walker can 
move only to the right, the probability of not making a 
step is 1/2, and therefore the correlation function 



C(t) = 2" 



(32) 



is purely exponential. 

In the general case < c < 1, the walker cannot leave 
the "cage" formed by two neighboring bridges. As al- 
ways, we consider the cage with sites x — 0, . . . ,£ — 1. 
For simplicity, let's assume again that the initial posi- 
tion is xq — 0. Rather than considering the walker in 
the cage (0, £ — 1) with a special behavior at x = and 
the absorbing boundary at x = £ one can analyze the or- 
dinary random walker in the extended cage (—£, £) with 
absorbing boundaries at x — £ and x = —£. The correla- 
tion function is merely the probability that this ordinary 
random walker will be at x = at time t and will remain 
inside the extended cage in intermediate times. This is a 
classic problem in probability theory whose solution is a 
cumbersome sum of expressions like (|30f) with alternat- 
ing (positive and negative) signs. Therefore we employ 
a continuum approximation which becomes asymptoti- 
cally exact when I 3> 1 (and accordingly c = t^ 1 <C 1). 
The solution is an infinite series of exponentially decay- 
ing terms. Keeping only the dominant term we obtain 



C{t) 



c exp 



7T 2 C 2 



t 



(33) 



as t — > oo. More precisely, the asymptotics (|33f) is valid 
when t ^> c~ 2 . In the regime 1 -C t <C c -2 , the domi- 
nant asymptotics is the same as in the c = case, that 
is, C(t) ~ i -1 / 2 . It is therefore understandable that a 
formula 



C(t)«(l + i) 



-1/2 



exp 



7T 2 C 2 



(34) 



fits well experimental data (and indeed it does 8]). Yet 
the true asymptotic behavior, Eq. (|33|l . is purely expo- 
nential without the power-law correction of Eq. (|34|l . 
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C. Modified burnt bridge model 

The precise definition of the walker dynamics at the 
boundary of the burnt bridge affects the results. We 
assumed that the walker at the boundary of the burnt 
bridge always moves to the the right. Recall, however, 
that in the stochastic version (p < 1) when the walker 
attempts to hop over the bridge from the left and the 
bridge burns, the walker actually remains at the same 
position. This suggests to modify the rule at the bound- 
ary of the burnt bridge — the walker either moves one 
step to the right or remains at the same position if it has 
tried the forbidden move across the burnt bridge. This 
defines the modified burnt bridge model. 

The calculation of v(c) and D{c) goes along the same 
lines as for the original burnt bridge model (Sec. IIII A|) 
and leads to the results presented in Sec. Inland displayed 
on Figs. 0and El 

For c = 1, the diffusion coefficient of Eq. (0} is the same 
D(l) — 1/8 in both the periodic and the random case. 
This particular result also follows from an independent 
calculation which we present here as it provides a good 
check of self-consistency. The key simplifying feature of 
the lattice fully covered with bridges is that the walker 
never hops to the left. The position Xt of the walker after 
t time steps satisfies 




probability 1/2 
probability 1/2 



from which 



and 



(xt+i) 



(Xt) 



<-*+!/ 



[Xt) 



(35) 



(36) 



(37) 



The variance a t = (x\) — (x t ) 2 satisfies a simple recur- 



rence 



0t+i = ot + j 



(38) 



which follows from Eqs. (|36|I - H37(I . The initial condition 
x a = implies (x ) = (Tq = 0. Solving (|3l)|) . subject 
to these initial values we obtain 



1 



(xt) 



1 



(39) 



The velocity and the diffusion coefficient can be read off 
the general relations (x t ) — vt and a t = 2Dt. Thus we 
recover the already known value v(l) = 1/2 and obtain 
the diffusion coefficient D(l) — 1/8 (which happens to 
be 4 times smaller than the bare diffusion coefficient). 



IV. THE STOCHASTIC BURNT BRIDGE 
MODEL 

Apart from randomness in hopping, the stochastic 
burnt bridge model has an additional stochastic element 
— crossing the bridge leads to burning with probability 
p while with probability 1 — p the bridge remains in- 
tact. To avoid the possibility of trapping we additionally 
assume that if the particle attempts to cross the bridge 
from the right and the bridge burns, the attempt is a fail- 
ure and the walker does not move. We have succeeded 
in computing v(c,p) in the situation when the bridges 
are equidistant. We again employ an approach involving 
auxiliary functions T(x) (see Appendix |A")> and L(x) (de- 
fined below). Perhaps, the entire problem can be treated 
by a direct approach discussed in Appendix [51 but that 
approach is more lengthy and we have only succeeded in 
computing the velocity for c = 1 that way. 

We must determine the average position of the first 
bridge that burns, and the average time of that event. 
The walker starts at x — 0, but it is again useful to con- 
sider a more general situation when the walker starts at 
an arbitrary position x. Denote by L(x) the average po- 
sition of the walker at the moment when the first bridge 
burns. The walker hops x — > x ± 1, and therefore 

L(x) = \\L{x-l)+L{x + l)\ (40) 

when x 7^ nl — 1, nl with n = 1,2,3,.... On the bound- 
aries of the bridges the governing equation (|4L)I) should 
be modified to account for possible burning events: 

L{ne _ 1)= L(ni-2) + (l-p)L(ni) +P ni ^ 
L(ni) = W + l) + (l-p)Hnt-l) + pnl ^ 

Equation l|4(J|) shows that L{x) is a linear function of x 
on each interval between the neighboring bridges, i.e., 



L(x) = A„, + (x - nl)B n 



(42) 



for nl < x < (n + 1)£ - 1. Plugging into 141a|) we 
obtain 

A„_! + £B n ^ = (1 - p)A n + pnl (43) 
Similarly, equation (|41bp reduces to 

A n = B n + (1 -p)[A n -x + {t-l)B n -x]+pnl (44) 
Using l|43l) . we get rid of £?'s in (|44|l and obtain 



A n -i -2gA n + A n+1 



l-p 

Here we used a shorthand notation 

p(2-p)(t-l) + 2 



P£ pni P+ } 2 - p)i (45) 



l-p 



fJ 



2(1 -P) 



(46) 
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The recurrence (|45|) admits a general solution 

A n = n£ + a + A+\ n + + A_ A™ (47) 

where A n = n£ + a with a = £/[p + (2 — p)f] is a par- 
ticular solution of the inhomogeneous equation 145(1 : the 
remaining contribution A + A™ + A_A™ with 



Ah 



9' 



1 



(48) 



is the general solution of the homogeneous part of l|45|) . 

If the walker is initially located far away from the ori- 
gin, x £, the first bridge would burn somewhere in its 
proximity, that is L(x) ~ x. This in conjunction with 
ijHJ lead to A n - n£ = 0(1) when n > 1. On the other 
hand, the general solution (|47|l grows exponentially since 
A+ > 1. This shows that the corresponding amplitude 
must vanish: A + = 0. Since L(x) is constant on the 
interval < x < I — 1, we have -Bo = 0, or [see (I43|) ] 



A = (1 - p)Ai + p£ 



By inserting Ao = a + A 
()49J) and solving for A 



we 



and A\ = £ 

jet 

^ — pa 
1-(1-P)A_ 



(49) 
A_A_ into 



(50) 



Return now to the situation when the walker starts at 
the origin. The average displacement of the walker after 
the first burning event is (x) = L(0) = Aq = a + A—, or 



p+(2-p)£ 



1 



(2-p)£ 



l-(l- P )X- 



(51) 



In the limiting cases p = 1 and £ — 1 we indeed recover 
(x) = £ and (|B9J| . respectively. 

The second part of the program is to determine the 
average time when the first burning occurs. Again we 
choose to investigate a more general quantity T(x). It 
satisfies Eq. I|A2[) when x ^ n£— 1, ni. On the boundaries 
of the bridges, the governing equations become 



T(n£- 1) 



T{n£-2) + (1 -p)T{n£) 



T(n£) = 



T(n£+1) + (l-p)T(n£-l) 



(52a) 
(52b) 



We seek a solution of l|A2(l . (|52a|) . il52b(l which is invari- 
ant under the transformation x <-» — x, and periodic in 
the large x limit. A solution to equation (|A2() is quadratic 
in x, viz. —x 2 + Yx + Z with arbitrary Y", Z. The same 
holds in our situation except that solutions in different 
intervals between the neighboring bridges differ. Thus 

TO) = ~{x-n£f + {x-n£)Y n + Z n (53) 

Plugging into l|52a|) - (|52bf) we obtain 

lY n -x + Z n -x = {l-p)Z n +£ 2 

Zn ~ l ~ Yn = Z n ^ + {£-l)Y n ^-(£-lf 
1 ~P 



Using the first equation, we exclude Y's from the second 
which turns into a recurrence 



7 7^7 , Q P + ( 2 ~ PY n 

Z n -\ - 2gZ n + Z n+ i + £ : = 



l-p 



whose general solution reads 



£ 
P 



Z,\l + Z_A™ 



(54) 



(55) 



The symmetry x <-> — x leads to Yq = 0, or 
Z = (1 - p)Z x + £ 2 



(56) 

The periodicity in the large x limit implies that Z n re- 
mains finite for large n. The exponentially growing part 
of the solution should therefore vanish, Z + = 0. Thus 
Z = Z_ +£/p and Z\ = Z_A_ +£/p. By inserting these 
relations into (|56|l we obtain 



£(£-!) 



l-(l-p)A_ 



(57) 



Thus the average time in the original problem is given by 
(t) = T(0) = Z , or 



<*> = p + r 



(58) 



(l-p)A- 

In the limiting cases p = 1 and £ = 1 we indeed recover 
(t) = f 2 and (t) = p- 1 [Eqs. (O and (JB2J)], respectively. 
Finally, the velocity is 

l-(l-p)A_ + (2-p)* 



P 



(59) 



(t) p+(2-p)€l-(l-p)A_+p(^-l) 

Using l|48[l . one can transform l|59|l into JSJ. 

In the initial state all bridges are intact, and a frac- 
tion of them remains intact as the walker moves along. 
The walker passes on average (x) / £ bridges per one burnt 
bridge. Hence the fraction / of bridges which forever re- 
mains intact approaches 

(x)/t-l 



I 



c)/l 



1 



(60) 



in the long time limit. Using Eq. (|51(l one recasts H60(l 
into Eq. ©. 

The calculation of the diffusion coefficient seems to be 
very challenging. For a system full of bridges (c = 1), 
however, the walk is somewhat analogous to the p = 1 
(and c < 1) case, and D might be possible to derive us- 
ing the approach presented in Sec. IIII Al The complete 
analysis appears to be very cumbersome, but if p — ► 
in addition to c = 1, the successive burnt bridges are 
(on average) separated by large gaps and therefore one 
can employ a continuous treatment. Following the steps 
described in Sec. IIII Al we obtained D = 1/4. This predic- 
tions agrees with simulations. Interestingly (see Fig. EJ, 
D is a non-monotonous function of p. 

Solving the stochastic burnt bridge model for randomly 
distributed bridges does not look possible in the realm of 
the above framework. Indeed, instead of working with 
ordinary deterministic recurrences like ()45|l one has to 
tackle stochastic recurrences (see Appendix [CJ . 
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FIG. 7: Numerical results for the diffusion coefficient for a 
system full of bridges (c = I) as function of p for several times. 
The arrow points to the analytic result D(p — > 0) = 1/4. 



V. DISCUSSION 

Our current understanding of the stochastic burnt 
bridge model is certainly incomplete — only the periodic 
case is somewhat tractable, albeit even in this situation 
we do not know how to compute various interesting quan- 
tities like the diffusion coefficient or the probability that 
in the final state two nearest burnt bridges are separated 
by k intact bridges. 

In many biological applications, molecular motors 
move along a homogeneous polymer filament (kinesin and 
myosin are classical examples jjj ) , while in other applica- 
tions the track is inhomogeneous (this particularly hap- 
pens when motors move along DNA). It would be inter- 
esting to study the burnt bridge model when in addition 
to the disorder related to location of the bridges there is 
the disorder associated with hopping rates. Earlier work 
on random walkers under the influence of a random force 
(El 01 ancl references therein) and recent work moti- 
vated by single-molecule experiments on motors moving 
along a disordered track 0, 0] may be useful in that 



regard. 
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APPENDIX A: CALCULATION OF [t] e AND [t\ 
FOR AN INTERVAL 

Let t be the first passage time, namely, the time it takes 
for the simple random walk in an interval [0, £] starting at 
site to reach site £ for the first time. Here we compute 



the first two moments, [t]e and [t 2 ]i, of this random vari- 
able. We will present an elementary approach that does 
not require the calculation of the complete first passage 
probability [ljj. (The calculations in Sec. lIVI can also be 
considered as a generalization of this method.) 

The process can be understood in terms of a random 
variable t(x), which is the time it takes to reach site £ 
if the walker starts at site x. As the walker from site x 
steps equally probably to either side 



t(x) 



t(x - 1) + 1 
t(x + l) + 1 



probability 1/2 
probability 1/2 



(Al) 



and for the average time T(x) = \t{x)]t we arrive at the 
recursion formula 

T(x) = ^[T(x-l)+T(x + l)] + l. (A2) 

This master equation holds for 1 < x < £ — 1 while for 
x = it should be replaced by 



T(0) = T(l) + 1 



(A3) 



since when the walker is at site 0, it always makes the step 
to the right. The recurrence (| A2|) — (| A3|> is supplemented 
by the boundary condition T{£) = 0. 

Equation (| A3|) can be re-written in the general form 
(|X2l) if T(-l) = T(l). Further one finds that Eq. i|Al)l 
holds for x = -1 if T(-2) = T(2), and generally Eq. l|Al|l 
applies for all \x\ < £ — 1. Hence we seek a solution in- 
variant under the transformation x *-* — x; the absorb- 
ing boundary conditions are T(±£) — 0. (Numerous ex- 
amples of analyzing equations like (| A2II with absorbing 
boundary conditions are described in 17].) The solution 
is very neat 



T(x) = (£ + x)(£-x) 



and in particular 



[t] t = t(o) = e 



(A4) 



(A5) 



For the derivation of the second moment, it is again 
convenient to consider T2{x) = [t 2 {x)]t which is the av- 
erage square time to reach site £ for the first time if the 
walker starts at position x. From Eq. IjAlfl one obtains 
the governing equation for 1 < x < £ — 1 



0= ^D 2 T 2 (x)+T(x- 



1) + T(x + 1) + 1 



(A6) 



where D 2 F(x) = F(x - 1) - 2F(x) + F(x + 1) is the 
shorthand notation for the discrete derivative of the sec- 
ond order. For x = we have 



T a (0) = r 2 (l) + 2T(1) + 1 



(A7) 



We can again seek a solution to Eq. (|A6|) satisfying 
the symmetry requirement x <-> — x and the absorbing 
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boundary conditions are T 2 (±£) = 0. Using Eq. (|A4|> we 
recast Eq. (|A6() into 

D 2 T 2 {x) = 4x(> + 1) -4.t + 2-4^ 2 , (A8) 

which yields to 

T 2 {x) = \x 2 [x 2 + 2-6£ 2 ] +[t 2 ] e 



with 



[t\ = if(5£ 2 - 2) 



(A9) 



The derivation of the moments for the modified burnt 
bridge model, where the hopping rule differs from the 
original model only from site 0, follows the same lines. 
The new rule affects only Eq. I|A3(1 which now becomes 
T(0) = |[T(Q) + T(l)] + 1. This equation can be recast 
in the general form (|A2|) if T(— 1) = T(0), and overall 
the symmetry T(x) = T(—x — 1) allows us to reduce 
the problem to solving i|A2|l subject to T{£) = and 
T(-£ - 1) = 0. The solution T(x) = (£ + 1 + x){£ - x) 
yields [t]t = T(0) = £{£+l). 

For the second moment the governing equation is given 
by Eq. for 1 < x < £ - 1, and for x — it is 

T 2 (0) = T2(0) + T2(1) + T(0) + T(l) + 1 . (A10) 

The boundary condition is T2(£) = 0. 

A solution of Eq. l|A8j) invariant under the transforma- 
tion x <-> —x — 1 and satisfying the absorbing boundary 
conditions T 2 (f) = T 2 (-£ - 1) = is 



D 2 T 2 (x) = Ax{x + 1) + 2 - 4(^ 2 + £) 
which is solved to yield 

T 2 (x) = -(x-l)x(x + l){x + 2) 



(All) 



with 



+ [1 - 2(£ 2 + £)] x[x + 1) + [t 2 ] e 



[t 2 ]e = -£(£ + l)[5£(£ + l)-l} 



(A12) 



APPENDIX B: DIRECT CALCULATION OF v(l,p) 

Here we present an alternative, direct calculation of the 
velocity for a lattice fully covered with bridges (c = 1). 
At each time step, the walker makes a move, so after 
t time steps all bridges remain intact with probability 
(1 — pY . Hence the first burning event would happen at 
time (t + 1) with probability 



B(t)= P (i- P y 



(Bl) 



and thus the average time till the first burning event is 
[t]=J2(t+l)p(l-p) t =p- 1 (B2) 



We also need the probability distribution P(x,t) of the 
position of the walker. As described earlier, we can con- 
sider the unconstrained random walk on the infinite line, 
and then "fold" it at the origin to give 



P(x,t) 



P (x,t) + P (-x,t) forx>0 



Po(0,i) 



for x = 



(B3) 



with Pq(x, t) being the probability distribution of the un- 
constrained walker. When the walker starts from the ori- 
gin at time t = 0, this probability is 

Pa(x,t)=< 2 "(4) i0Tt + XGVen (B4) 
for t + x odd 

The probability that a bridge burns at time (t + 1) when 
the walker hops from site x is P(x,t)B(t), and the total 
probability is obtained after summing over all t: 



%(x) = J2P(x,t)B(t) . 



(B5) 



t=o 



If the walker is hopping to the right when the burn- 
ing occurs, the move is completed; otherwise the walker 
remains in its position. Both of these alternatives occur 
equiprobably when x > 0, while when x — the walker 
surely hops to the right. The average final position of the 
walker is therefore 



23(0)+^ 



x=l 



Z+i)'fM.r) |B(i) 



Using l|B3|) . l|B5|) . and the identity J2 x >q' B ( x ) = 1 we 
transform l|B6() into 



2 2 



t=o 



x — — oo t— 



The first sum reduces to 

oo 



k=0 



(B7) 



2f- V k J 2V1 - 4a 2 
where a = (1 — p)/2. Next we re- write the second sum as 

oo t 

X>(W) . v(t)= J2 M-Po(m) 

f=0 x=-t 

and simplify V(t) by separately considering even and odd 
times: 



2 -2fc+l ^ | m | / 



2fe \ 
k+mJ 



for t = 2fc 



2- 2fe E H O for* = 2fc + l 

m— — k 
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Evaluating the sums we obtain 

2k 

| 2-^(fc + l)( ^ 

V(t) 



2- 2k {k + l) 



\k + 1 
2k + 1 



k + 1 



for i = 2k 
for i = 2fc + 1 



Putting this into J2t>o we nn d that the sum is 

equal to 



2px; 



fc=0 

2p 



2a 2 + a 
(1 - 4a 2 ) 3 / 2 



(B8) 



Combining (|B7J| and (|BS|) we obtain the average displace- 
ment 



1 1 2-p 



(B9) 



2 2 V p 

and therefore w = (x)/[t] = p(x) is indeed given by © 



APPENDIX C: THE STOCHASTIC BURNT 
BRIDGE MODEL IN THE CASE OF RANDOMLY 
POSITIONED BRIDGES 



The formalism detailed in Sect. II VI for the periodic lo- 
cation of bridges to the situation formally applies to the 
situation when bridges are arbitrarily distributed. Let 
(£% — l,£i), (£i +£2 — +£2), etc. be bridge locations. 
Away from bridges the governing equation (|4(JI) is valid 
while on the boundaries the modified equations are al- 
most identical to (|41a|l and l|41b(l . the only exception is 
that n£ should be replaced by L n = £\ + . . . + £ n . The 
average position of the walker L(x) is again a linear func- 
tion of x on each interval between neighboring bridges; 



for L n < x < L n + £ n+ i - 1 

L(x) — A n + (x - L n )B r , 
The analogs of Eqs. !)33 )) -l|13 jl are 



(CI) 



B n -i = [-A,-! + {l-p)A n +pL n ]/£ n 

A n = B n + (1 -p)[A n ^l + (£ n - l)B n -l] + pL n 

Using the first equation we exclude £Ts from the second 
and thereby recast it into a recurrence 



A 



n+l 



A r . 



An 



£ n 



+1 



l-p 
+ p + pL 



' , n n 1 (!-p) 2 

p{2 - p) + £— + -T^ 



2-p 



1 



l-p 



£m 



f 



(C2) 



In the interesting case when fs are independent iden- 
tically distributed random variables, one must solve the 
stochastic inhomogeneous recurrence l|C2ll . Even a ho- 
mogeneous version of Eq. i|C2(l is analytically intractable. 
The additional challenging feature of the inhomogeneous 
recurrence (|C2|) is infinite memory manifested by factor 
L n ; as a result, it is not clear how to find a particular 
solution of Eq. I|C2(1 which is required if one wants to 
reduce the problem to solving a homogeneous version of 
Eq. 

The case of weak disorder is probably exceptional, e.g., 
it should be possible to compute the growth rates X±. 
One can, however, avoid such a lengthy analysis by not- 
ing that in the present context the condition of weak 
disorder implies that bridges are located almost period- 
ically and their concentration is small (c <C 1). Assum- 
ing additionally that the bridge burning probability is 
not anomalously small, so that c <SC p, an argument pre- 
sented after Eq. (J7J) shows that in the leading order the 
burnt bridge model must be recovered. Thus v ~ c and 
D w (l-c 2 )/3. 
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